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Abstract 

We consider a system of nonlinear partial differential equations with stochastic 
dynamical boundary conditions that arises in models of neurophysiology for the 
diffusion of electrical potentials through a finite network of neurons. Motivated by 
the discussion in the biological literature, we impose a general diffusion equation 
on each edge through a generalized version of the FitzHugh-Nagumo model, while 
the noise acting on the boundary is described by a generalized stochastic Kirchhoff 
law on the nodes. In the abstract framework of matrix operators theory, we rewrite 
this stochastic boundary value problem as a stochastic evolution equation in infinite 
dimensions with a power-type nonlinearity, driven by an additive Levy noise. We 
prove global well-posedness in the mild sense for such stochastic partial differential 
equation by monotonicity methods. 

1 Introduction 

In this paper we study a system of nonlinear diffusion equations on a finite network 
in the presence of an impulsive noise acting on the nodes of the system. We allow a 
rather general nonlinear drift term of polynomial type, including functions of FitzHugh- 
Nagumo type (i.e. f{u) = —u{u — 1) (n — a)) arising in various models of neurophysiology 
(see e.g. the monograph [19] for more details). 

Electric signaling by neurons has been studied since the 50s, starting with the now 
classical Hodgkin-Huxley model [16] for the diffusion of the transmembrane electrical po- 
tential in a neuronal cell. This model consists of a system of four equations describing the 
diffusion of the electrical potential and the behaviour of various ion channels. Successive 
simplifications of the model, trying to capture the key phenomena of the Hodgkin-Huxley 
model, lead to the reduced FitzHugh-Nagumo equation, which is a scalar equation with 
three stable states (see e.g. [27]). 
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Among other papers dealing with the case of a whole neuronal network (usually 
modeled as a graph with ra edges and n nodes), which is intended to be a simplified 
model for a large region of the brain, let us mention a series of recent papers by Mugnolo 
et al. [211|2l], where the well-posedness of the isolated system is studied. 

Note that, for a diffusion on a network, other conditions must be imposed in order to 
define the behaviour at the nodes. We impose a continuity condition, that is, given any 
node in the network, the electrical potentials of all its incident edges are equal. Each 
node represents an active soma, and in this part of the cell the potential evolves following 
a generalized Kirchhoff condition that we model with dynamical boundary conditions 
for the internal dynamics. 

Since the classical work of Walsh [28], stochastic partial differential equations have 
been an important modeling tool in neurophysiology, where a random forcing is intro- 
duced to model external perturbations acting on the system. In our neuronal network, 
we model the electrical activity of background neurons with a stochastic input of impul- 
sive type, to take into account the stream of excitatory and inhibitory action potentials 
coming from the neighbors of the network. The need to use models based on impul- 
sive noise was already pointed out in several papers by Kallianpur and coauthors - see 
e.g. [17[ I18| . On the other hand, from a mathematical point of view, the addition of 
a Brownian noise term does not affect the difficulty of the problem. In fact, in section 
[2] below, a Wiener noise could be added taking (7 7^ 0, introducing an extra term that 
does not modify the estimates obtained in section [Sj which are the basis for the principal 
results of this paper. Let us also recall that the existence and uniqueness of solutions 
to reaction-diffusion equations with additive Brownian noise is well known - see e.g. 

Following the approach of we use the abstract setting of stochastic PDEs by 
semigroup techniques (see e.g. [lOlIII]) to prove existence and uniqueness of solutions 
to the system of stochastic equations on a network. In particular, the specific stochastic 
dynamics is rewritten in terms of a stochastic evolution equation driven by an additive 
Levy noise on a certain class of Hilbert spaces. 

The rest of the paper is organized as follows: in section [2] we introduce the problem 
and we motivate our assumptions in connection with the applications to neuronal net- 
works. Then we provide a suitable abstract setting and we prove, following j2^, that 
the linear operator appearing as leading drift term in the stochastic PDE generates an 
analytic semigroup of contractions. Section [3] contains our main results. First we prove 
existence and uniqueness of mild solution for the problem under Lipschitz conditions 
on the nonlinear drift term. This result (essentially already known) is used to obtain 
existence and uniqueness in the mild sense for the SPDE with a locally Lipschitz drift 
of FitzHugh-Nagumo type by monotonicity techniques. 

2 Setting of the problem 

Let us begin introducing some notation used throughout the paper. We shall denote 
by ^ and respectively, weak and weak* convergence of functions. All stochastic 
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elements are defined on a (fixed) filtered probability space ($7,^, JF^jP) satisfying the 
usual hypotheses. Given a Banach space E, we shall denote by U'{E) the space of 
-E- valued random variables with finite p-th moment. 

The network is identified with the underlying graph G, described by a set of n 
vertices vi, . . . , v„ and m oriented edges ei, . . . , which we assume to be normalized, 
i.e., ej = [0,1]. The graph is described by the incidence matrix ^ = — , where 
= i^ti)nxm and = (0j^)„xm are given by 

Vi = ej(l) ^+ ^ A' Vi = ej(0) 
otherwise lO, otherwise. 

The degree of a vertex is the number of edges entering or leaving the node. We denote 

r(vi) = {j £ {1,... ,m} : ej(0) = Vj or 6^(1) = v^} 

hence the degree of the vertex Vj is the cardinality |r(vj)|. 

The electrical potential in the network shall be denoted by u{t, x) where u G (^'^(0, 1))" 
is the vector {ui{t, x), . . . , Um{t, x)) and Uj{t, •) is the electrical potential on the edge ej. 
We impose a general diffusion equation on every edge 

d d f d \ 

= ^ lCj{x)—Uj{t,x) j + fj{Uj{t,x)), (1) 

for all (t,x) E x (0,1) and all j = l,...,m. The generality of the above diffusion 
is motivated by the discussion in the biological literature, see for example [T9], who 
remark, in discussing some concrete biological models, that the basic cable properties is 
not constant throughout the dendritic tree. The above equation shall be endowed with 
suitable boundary and initial conditions. Initial conditions are given for simplicity at 
time t = of the form 

uj{0,x) =Ujo{x) eC{[0,l]), j = l,...,m. (2) 

Since we are dealing with a diffusion in a network, we require first a continuity assump- 
tion on every node 

Pi{t) := Uj(t,Vi) = Uk{t,Vi), t > 0, i, e r(vi), i = 1, ...,n (3) 
and a stochastic generalized Kirchhoff law in the nodes 

^Pi{t) = -biPiit) + ^ (j)ijHjCj(yi)^Uj{t,Vi) + ai-^L{t,yi), (4) 
i6r(v,) 

for all t > and i = 1, . . . ,n. Observe that the plus sign in front of the Kirchhoff term 
in the above condition is consistent with a model of purely excitatory node conditions, 
i.e. a model of a neuronal tissue where all synapses depolarize the postsynaptic cell. 
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Postsynaptic potentials can have graded amplitudes modeled by the constants fij > 
for all j = 1, m. 

Finally, L{t,Vi), i = l,...,n, represent the stochastic perturbation acting on each 
node, due to the external surrounding, and -^L(t,Vi) is the formal time derivative of 
the process L, which takes a meaning only in integral sense. Biological motivations 
lead us to model this term by a Levy process. In fact, the evolution of the electrical 
potential on the molecular membrane can be perturbed by different types of random 
terms, each modeling the influence, at different time scale, of the surrounding medium. 
On a fast time scale, vesicles of neurotransmitters released by external neurons cause 
electrical impulses which arrive randomly at the soma causing a sudden change in the 
membrane voltage potential of an amount, either positive or negative, depending on the 
composition of the vesicle and possibly even on the state of the neuron. We model this 
behaviour perturbing the equation by an additive n-dimensional impulsive noise of the 
form 



See Hvpothesis 12.21 below for a complete description of the process and [TS] for a related 
model. 

Although many of the above reasonings remain true also when considering the dif- 
fusion process on the fibers, we shall not pursue such generality and assume that the 
random perturbation acts only on the boundary of the system, i.e. on the nodes of the 
network. 

Let us state the main assumptions on the data of the problem. 
Hypothesis 2.1. 

1. In |2P, we assume that Cj(-) belongs to C^([0, 1]), for j = 1, . . . ,m and Cj{x) > 
for every x € [0, 1]. 

2. There exists constants r/ € M, c > and s > 1 such that, for j = l,...,m, 
the functions fj{u) satisfy fj{u) + rju is continuous and decreasing, and \fjiu)\ < 
c(l + 

3. In Q), we assume that bi > for every i = l,...,n and at least one of the 
coefficients hi is strictly positive. 

4- {fJ-j}j=i,...,m CLnd {(Tj}i=i,...,n o'^c real positive numbers. 

Given a Hilbert space 7i, let us define the space L%{n X [0, T] ; of adapted processes 
y : [0, T] ^ endowed with the natural norm 



\ Jo J 

We shall consider a Levy process {L{t), t > 0} with values in (M",;B(M")), i.e., a 
stochastically continuous, adapted process starting almost surely from 0, with stationary 




(5) 
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independent increments and cadlag trajectories. By the classical Levy-Ito decomposition 
theorem, one has 

L{t) = mt + qWt + x[N{t,dx) -tu{dx)]+ xN{t,dx), t>0 (6) 

J\x\<l J\x\>l 

where m £ M", q € MnxnO^) is a symmetric, positive defined matrix, {Wt, t > 0} is an 
n-dimensional centered Brownian motion, N{t, dx) is a Poisson measure and the Levy 
measure i/(dx) is cj-finite on M" \ {0} and such that J min(l, x^)z^(dx) < oo. We denote 
by N{dt,dx) := N{dt,dx) — dtz^(dx) the compensated Poisson measure. 

Hypothesis 2.2. We suppose that the measure v has finite second order moment, i.e. 

|xpi/(dx) < oo. (7) 



Condition ([7]) implies that the generalized compound Poisson process J^^i x N{t, dx) 
has finite moments of first and second order. Then, with no loss of generality, we assume 
that 

/ xi^{dx) = 0. (8) 

^|x|>l 

We also assume throughout that the Levy process is a pure jump process, i.e. m = 
and q = 0, which leads to the representation ([5]) in view of assumptions ([7|) and ([8]). 

2.1 Well-posedness of the linear deterministic problem 

We consider the product space IH = (L^(0, 1))™. A vector -u € EI is a collection of 
functions {uj{x), x G [0,1], j = l,...,m} which represents the electrical potential 
inside the network. 

Remark 2.3. For any real number s > we define the Sobolev spaces 

= (F^(o,i)r, 

where -ff*(0, 1) is the fractional Sobolev space defined for instance in [22\. In particular 
we have that C (C[0, 1])"^. Hence we are allowed to define the boundary evaluation 
operator n : H-'^ ^ M" defined by 



Uu 



, where Pi = li(vi) = tifc(vi) for /c G r(vi), i = 1, n. 

\Pn^ 

On the space M we introduce the linear operator (A, D{A)) defined by 



D{A) = {u £lf \3p eR" such that Uu = p} 
— (cix) — 

dx ^ dx' J J j=l,...,m 
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As discussed in [23], the diffusion operator A on a network, endowed with active 
nodes, fits the abstract mathematical theory of parabolic equations with dynamic bound- 
ary conditions, and in particular it can be discussed in an efficient way by means of 
sesquilinear forms. 

Notice that no other condition except continuity on the nodes is imposed on the 
elements of D{A). This is often stated by saying that the domain is maximal. 

The so called feedback operator, denoted by C, is a linear operator from D(A) to 
defined as 

/ 

Cu= [ (l)ijHjCj{vi)—Uj{t,\/i) 

Vier(v,) 



i=l,...,n 



On the vector space M" we also define the diagonal matrix 



B 



\ 



With the above notation, problem ([T])-([l]) can be written as an abstract Cauchy 
problem on the product space 7^ = H x endowed with the natural inner product 

(X, Y)-u = {u, v)u + {p, q)R" , where X,Y eH and X = [ ] , Y 



We introduce the matrix operator A on the space 7Y, given in the form 



A 



A 
C B 



with domain 

D{A) = {X = {u,p) £n : ue D{A),Uj{yi) = pi for every j G r(v,)}. 
Then the linear deterministic part of problem becomes 

ixit) = AXit) 
X(0) = xo 

where xq = (uj(0, x))j=i_...^m € C([0, 1])™ is the vector of initial conditions. This problem 
is well posed, as the following result shows. 

Proposition 2.4. Under Hypotheses \2.1^ 1\ and \2.1^ the operator {A,D{A)) is self- 
adjoint, dissipative and has compact resolvent. In particular, it generates a Cq analytic 
semigroup of contractions. 
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Proof. For the sake of completeness, we provide a sketch of the proof following [24J. The 
idea is simply to associate the operator {A,D{A)) with a suitable form a{X,Y) having 
dense domain V CTC. 

The space V is defined as 

V=|x= Q \ ue{HH0,l)r, Ukivi) = Pi for i = l,...,n,keT{yi)^ 

and the form o is defined as 



3=1-'^ i=\ 



The form a is clearly positive and symmetric; furthermore it is closed and continuous. 
Then a little computation shows that the operator associated with a is -D(^)) defined 
above. Classical results in Dirichlet forms theory, see for instance [25], lead to the desired 
result. □ 

The assumption that 6; > for some / is a dissipativity condition on A. In particular 
it implies the following result (for a proof see |24j ) . 

Proposition 2.5. Under Hypotheses \2. il fTl and \2.1t ^ the operator A is invertible and 
the semigroup {T{t), t > 0} generated by A is exponentially bounded, with growth bound 
given by the strictly negative spectral bound of the operator A. 



3 The stochastic Cauchy problem 

We can now solve the system of stochastic differential equations ([T])- dH) . The functions 
fj{u) which appear in ([1]) are assumed to have a polynomial growth. We remark that 
the classical FitzHugh-Nagumo problem requires 

fj{u) = u{u-l){aj-u) i = l,...,m 

for some aj € (0, 1), and satisfies Hvpothesis 12.1121 with 

+ 1) 



We set 



F{u) = {f,{n,))^^^_^ and ^(X) = ( ^^^^^^ for X = (^^ , (10) 



and we write our problem in abstract form 



j dX{t) = [AX{t) - T{X{t))] dt + S dC{t) 
[ X{0) = xo, 



ill) 
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where S is the matrix defined by 



[O crj [O diag(cri, . . . ,cj„; 
and C{t) is the natural embedding in Ti of the n-dimensional Levy process L{t), i.e. 

Remark 3.1. Note that in general is only defined on its domain D{J^), which is strictly 
smaller than 7i. 



Let us recall the definition of mild solution for the stochastic Cauchy problem (jllj) . 

Definition 3.2. An 7i-valued predictable process X{t), t € [0, T], is said to be a mild 
solution of ( fTTj) if 

[ \J^{Xis))\ds <+oo (12) 

(13) 



and 



X{t) = T{t)xQ- I T{t- s)T{X{s))ds+ [ T{t-s)^dC{s) 
Jo Jo 

P-a.s. for all t € [0,T], where T{t) is the semigroup generated by A. 

Condition p2p implies that the first integral in (|13p is well defined. The second 
integral, which we shall refer to as stochastic convolution, is well defined as will be 
shown in the following subsection. 

3.1 The stochastic convolution process 

In our case the stochastic convolution can be written as 



The definition of stochastic integral with respect to a compensated Poisson measure 
has been discussed by many authors, see for instance [H EJ [3l [U [HI [15]. Here we 
limit ourselves to briefly recalling some conditions for the existence of such integrals. 
In particular, in this paper we only integrate deterministic functions, such as T(-)E, 
taking values in (a subspace of) L{Ti.), the space of linear operators from H to TC. In 
order to deflne the stochastic integral of this class of processes with respect to the Levy 
martingale- valued measure 

M{t,B) = [ xN{t,dx), (14) 
Jb 
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one requires that the mapping T(-)S : [0,T] x M" 9 {t,x) i-^ T{t){0,ax) belongs to the 
space L2((0,r) x B; (M(dt,dx))) for every B G i3(M"), i.e. that 



T 




Tis) 







i^(dx) ds < oo. 



(15) 



H 



Thanks to ([7]), one has 



JB 



T{s) 



< a 





ax 

T 



^{dx) ds 



n 



\ns)\ 







L(n) 



ds 



jxp z^(dx) I < oo, 



B 



thus the stochastic convolution Z{t) is well defined for all t G [0, T]. 

We shall now prove a regularity property (in space) of the stochastic convolution. 
Below we will also see that the stochastic convolution has cadlag paths. 

Let us define the product spaces £ := (C[0, 1])™ x M" and C:f{[0,T];L^{n;£)), the 
space of £'-valued, adapted mean square continuous processes Y on the time interval 
[0, T] such that 

\Y\l^ := sup E\Y{t)\l <oo. 

t£[0,T] 

Lemma 3.3. For all t G [0, T], the stochastic convolution {Z{t), t G [0,T]} belongs to 
the space Cjr{[0,T]; L'^{Q;£)). In particular, Z{t) is predictable. 

Proof. Let us recall that the (unbounded) matrix operator ^ on 7-^ is defined by 

^_(dl 
"^"V-^. B 

with domain D{A) = {X = {u,p) € Ti. : u & D{A),ui(\/i) = Pi for every I G r(vj)}, and, 
by proposition 12.41 it generates a Co-analytic semigroup of contractions on H.. 

Let us introduce the interpolation spaces Tie = {'H,D{A))e.2 foi" ^ € (Ojl)- By 
classical interpolation theory (see e.g. [2^) it results that, for 6 < 1/4, = x M" 
while for 9 > 1/4 the definition of Tig involves boundary conditions, that is 



Hg 



e H"^^ X : Uu = p 



Therefore, one has {0,ax) G Tig for 9 < 1/4. Furthermore, for 9 > 1/2, one also has 



Hg C 



C (C[0, 1])™ X M" by Sobolev embedding theorem. Moreover, for all 



X eHg and 6* + 7 G (0, 1), it holds 

\T{t)x\g+^ < t-^\x\ge^^\ 
where io^ is the spectral bound of the operator A. 
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Let (9,7 be real numbers such that 9 G (0, 1/4), 7 G (0, 1/2) and 6* + 7 G (1/2, 1). 
Then for all t G [0, T] 



T(t - s) 



ax 



:"-a.s. 



?+7 



The right hand side of the above inequality is well defined if and only if 



E 



T{s) 





ax 



N{dx,ds) 



H7 



T 



T{s) 





ax 



i/(dx)ds < 00, 



H7 



where the identity follows by the classical isometry for Poisson integrals. On the other 
hand, one has 



T 



T{s) 





ax 



2 i-T 
i/(dx)ds < 

e+7 







-27 





ax 



22'^^V(da;)ds 



Jo JR" 



using 7 G (0,1/2) and assumption ([7]). So Z(t) G He+^y for + 7 > 1/2 and then 
Z(t) G (C[0, l])'" X M" = It remains to prove that Z{t) is mean square continuous as 
^-valued process. For 0<s<t<Twe can write 

2 



E\Z{t)-Zis)\i = E 
< 2E 
+2E 



[ T(t-r)Sd£(r) - / T{s-r)mC{r) 
Jo Jo 



[r{t-r)-T{s-r)] 




ax 



iV(dx,dr) 



T{t - r) 



ax 



N{dx, dr) 



[T{t - r) - T{s - r)] 





ax 



i/(dx)dr 



+2 



T(t - r) 





ax 



i^(dx)dr — > 



by the strong continuity of the semigroup T{t). Since the stochastic convolution Z{t) is 
adapted and mean square continuous, it is predictable. □ 



3.2 Existence and uniqueness in the Lipschitz case 

We consider as a preliminary step the case of Lipschitz continuous nonlinear term and 
we prove existence and uniqueness of solutions in the space Cjf of adapted mean square 
continuous processes taking values in 7i. We would like to mention that this result is 
included only for the sake of completeness and for the simplicity of its proof (which is 
essentially based only on the isometry defining the stochastic integral). In fact, a much 
more general existence and uniqueness result was proved by Kotelenez in [20j . 
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Theorem 3.4. Assume that Hypothesis \2.2\ holds, and let xq be an J^Q-measurable 7i- 
valued random variable such that E|xoP < oo. Let G : TC ^ Ti. be a function satisfying 
Lipschitz and linear growth conditions: 

|G(x)| <co(l + |x|), \G{x) - G{y)\ < co\x - y\, x,y dU. (16) 

for some constant cq > 0. Then there exists a unique mild solution X G C([0, T]; L^(i7, Ti)) 
to equation Ml\) with —T replaced by G. Moreover, the solution map xq X(t) is Lip- 
schitz continuous. 

Proof. We follow the semigroup approach of |1H Theorem 7.4] where the case of Wiener 
noise is treated. We emphasize only the main differences in the proof. 

The uniqueness of solutions reduces to a simple application of Gronwall's inequal- 
ity. To prove existence we use the classical Banach's fixed point theorem in the space 
Gjr{[0,T];L'^{n;n)). Let JC be the mapping 

/C(y)(t) = T{t)xo + /* T{t - s)G{Y{s)) ds + Z{t) 
Jo 

where Y G Cjr([0, T]; 7^)) and Z{t) is the stochastic convolution. Z{-) and T(-)xo 

belong to Cyr{[0, T];L^{n; H)) respectively in view of Lemma [3.3l and the assumption on 
Xq. Moreover, setting 

/Ci(y)(t) = / T{t-s)G{Y{s))ds, 
Jo 

it is sufficient to note that 

\JCi{Y)\l^<{Tcof{l + \Y\ 

by the linear growth of G and the contractivity of Then we obtain that fC maps 

the space Cjr{[0,T];L'^{Q.;n)) to itself. Furthermore, using the Lipschitz continuity of 
G, it follows that for arbitrary processes li and Y2 in Cjf([0, T]; 7^)) we have 

|/c(yi) - JC{Y2)\l^ = i/Ci(yi) - /Ci(y2)i^^ < {c^Tf\Y^ - Y^fc^. 

If we choose an interval [0,7"] such that T < Cq^, it follows that the mapping K, has 
a unique fixed point X G Cjr([0, T]; 7^)). The extension to an arbitrary interval 

[0, T\ follows by patching together the solutions in successive time intervals of length T . 

The Lipschitz continuity of the solution map xq 1-^ ^ is again a consequence of 
Banach's fixed point theorem, and the proof is exactly as in the case of Wiener noise. 

It remains to prove the mean square continuity of X. Observe that T(-)xo is a deter- 
ministic continuous function and it follows, again from Lemma 13.31 that the stochastic 
convolution Z{t,^ is mean square continuous. Hence it is sufficient to note that the same 
holds for the term ^^T{t — s)G{X{s)) ds, that is P-a.s. a continuous Bochner integral 
and then continuous as the composition of continuous functions on [0,T]. □ 

Remark 3.5. By standard stopping time arguments one can actually show existence and 
uniqueness of a mild solution assuming only that 2:0 is .Fo-measurable. 
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In order to prove that the solution constructed above has cadlag paths, unfortunately 
one cannot adapt the factorization technique developed for Wiener integrals (see e.g. 
[llj). However, the cadlag property of the solution was proved by Kotelenez [20], under 
the assumption that A is dissipative. Therefore, thanks to proposition 12.41 the solution 
constructed above has cadlag paths. One could also obtain this property proving the 
following a priori estimate, which might be interesting in its own right. 

Theorem 3.6. Under the assumptions of theorem \3.4\ the unique mild solution of prob- 
lem (jlip verifies 

E sup \X{t)\j^ < oo. 
te[o,r] 

Proof. Let us consider the Ito formula for the function | • |^, applied to the process 
X. Although our computations are only formal, they can be justified using a classical 
approximation argument. We obtain 

d|X(t)|2, = 2{Xit-),dXit))n + d[X]it). 

By the dissipativity of the operator A and the Lipschitz continuity of G, we obtain 

{X{t-),dX{t))n = {AX{t),X{t))ndt+{G{X{t)),X{t))ndt+{X{t-),mC{t))n 
<co\X{t)\l^ + {X{t-),mC{t))n. 

Therefore 

\X{t)\l^<\xo\lc + 2co f \X{s)\lcds + 2 [\x{s-),mC{s))n+ [' m^Cjis) 
Jo Jo Jo 



and 



E sup \X{t) <E|xo|>^ + 2coTE sup \X{t)\j^ 

t<T t<T 



+ 2E sup 

t<T 

where we have used the relation 



/ {X{s-),mC{s))n +T [ m^\x\^u{dx), (17) 
Jo Jr" 



Esup[X](t) <E / lSpd[£](t) =E /" \i:\^d{C){t) = T [ 

t<T Jo Jo Jk' 



0^ ' 

X 



v{dx). 



By the Burkholder-Davis-Gundy inequality applied to Mt = Jq{X{s—), Sd>C(s))7^, there 
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exists a constant ci such that 



E sup 

t<T 



{X{s-),m£{s))n 



< ciE 



{X{s-),mC{s))n 



(T) 



1/2 



1/2 



< ciE (sup|X(t)||^ / |Spd[/:](s) 

yt<T Jo 

< ci (eEsup|X(t)|^ + -J-E / |Epd[>C](s) I 

y t<T 4e Jq J 

= cieEsup|X(t)|f^ + ^ [ |S|2|xpiy(dx), (18) 



where we have used the elementary inequahty ab < ea^ + Then by ([T7j) and ([18 

we have 

Esup|X(t)|f^ < E\xo\n + 2coTEsup\X{t)\li^ + 2cieEsnp\X{t)\li 

t<T t<T t<T 



+ (^ + l)r / m'^\x\^iy{dx), 



hence 



where 



Esup|X(t)|^ <iv[E|xo|?^ + rfl + ;^ 



< +00, 



N = Nico,ci,T,e) 



1 



1 - 2coT - 2cie 

Choosing e > and T > such that < 1, one obtains the claim for a small time 
interval. The extension to arbitrary time interval follows by classical extension argu- 
ments. □ 

3.3 FitzHugh-Nagumo type nonlinearity 

Let us now consider the general case of a nonlinear quasi-dissipative drift term J^. 
Theorem 3.7. Let T : D{T) C 7i ^ TC be defined as in ilO\}. Then the equation 



j dX{t) = [AX{t) - J^{X{t))]dt + ^dC{t), t G [0,T], 
\ ^(0) = xo 

admits a unique mild solution, denoted by X{t,XQ), which satisfies the estimate 

E\X{t,x) - X{t,y)\'^ < e2^*E|x-yp. 

for all X, y G H. 



(19) 
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Proof. As observed in section [3] above, there exists 77 > such that F + rjl is accretive. 
By a standard argument one can reduce to the case of ?] = (see e.g. 0]), which we 
shall assume from now on, without loss of generality. Let us set, for A > 0, Fx{u) = 
F((l + XF)~^{u)) (Yosida regularization). J^x is then defined in the obvious way. 

Let Qy = —Ay + J'{y)- Then Q is maximal monotone on Ti. In fact, since A is 
self-adjoint, setting 

^\A^'M\ ueD{A^'^) 
+00, otherwise, 



one has A = dip. Let us also set F = dg, where g : — > M is a convex function, the 
construction of which is straightforward. Well-known results on convex integrals (see 
e.g. [H sec. 2.2]) imply that F on H is equivalently defined as F = dig, where 

[ / 9iuix))dx, if giu) € LH[0,ir), 
Ig{u) = I J[o,i]- 

[^+00, otherwise. 

Let us recall that 

Since D{J^) n D{A) is not empty, Q is maximal monotone if ip{{I + A^) ^(u)) < ip{u) 
(see e.g. [6l Thm. 9]), which is verified by a direct (but tedious) calculation using the 
explicit form of A, since (I + Xfj)~^ is a contraction on M for each j = 1, . . . , m. 
Let us consider the regularized equation 

dXx{t) + gxXx{t)dt = j:dCit). 

Appealing to Ito's formula for the square of the norm one obtains 

\Xx{t)f + 2 [\gxXx{s),Xx{s))ds = \xo\^ + 2 [\xx(,s-),J:dC{s)) + [XxM 
Jo Jo 

for all t E [0, T]. Taking expectation on both sides yields 

E\Xx{t)\^ +2E [\gxXx{s),Xx{s))ds = \xo\^ + t [ m^\z\^uidz), (20) 

Jo JR" 
where we have used the identity 

E[Xx]it)=E\ [ i:dC{s)\{t) = t [ |Sp|z[2z^(dz). 
L Jo JR" 

Since by 1^ we have that {Xx} is a bounded subset of L°°{[0,T],L'^{n)), and L^C^) 
is separable, Banach-Alaoglu's theorem implies that 

Xx^X inL^{[0,T],L\n)), 
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on a subsequence still denoted by A. Thanks to the assumptions on /-,-, one can easily 
prove that {F{u),u) > c|u|^+"'^ for some c > and p > 1, hence ()20p also gives 

E r \Xx{s)\lXlds<C, 
Jo 

which implies that 

Xx^X mLP+^{nx[0,T]x D,FxdtxdC), (21) 
where D = [0, l]™ x M"'. Furthermore, ([20]) and ([2T]) also imply 

GxXx^V mL^{nx[0,T]xD,¥xdtxdO. 

The above convergences immediately imply that X and rj are predictable, then in order 
to complete the proof of existence, we have to show that r]{uj,t,^) = Q{X{uj,t,^)), 
¥ X dt X dS^-a.e.. For this it is enough to show that 

limsupE /" {gxXx{s),Xxis))ds <E [ {r]{s), X{s)) ds. 

A-»0 Jo Jo 

Using again Ito's formula we get 

E\X{T)\^ + 2e[ {r]{s),X{s)) ds = IxqI"^ + T [ |Sp (22) 
Jo Jr" 

However, (j2T]) implies that 

liminf E|XA(r)P > mX(T)\^ 

(see e.g. [TJ Prop. 3.5]), from which the claim follows comparing ()20p and (I22p . 

The Lipschitz dependence on the initial datum as well as (as a consequence) unique- 
ness of the solution is proved by observing that X{t,x) — X(t,y) satisfies P-a.s. the 
deterministic equation 

j^{X{t, x) - X{t, y)) = A{X{t, x) - X{t, y)) - T{X{t, x)) + T{X{t, y)), 

hence 

i^|X(t,x)-X(t,y)|2 = {A{X{t,x)-X{t,y)),X{t,x)-X{t,y)) 

-{HX{t, x) - J^{X{t, y)),X{t, x) - X{t, y)) 
< r]\X{t,x)-X{t,y)\^, 

where X{-,x) stands for the mild solution with initial datum x. By Gronwall's lemma 
we have 

E|X(t,x) -X(t,y)|2 < e2^*E|x-yp, 
which concludes the proof of the theorem. □ 
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Remark 3.8. An alternative method to solve stochastic evolution equations with a dis- 
sipative nonlinear drift term is developed in |1H I12j . for the case of Wiener noise, and 
in the recent book [26] for the case of Levy noise. This approach consists essentially in 
the reduction of the stochastic PDE to a deterministic PDE with random coefficients, 
by "subtracting the stochastic convolution". To carry out this plan one has to find a 
reflexive Banach space V, continuously embedded in 7i, which is large enough to contain 
the paths of the stochastic convolution, and at the same time not too large so that it 
is contained in the domain of the nonlinearity J-. In particular, in the case of equation 
([19]), theorem 10.14 in [26] yields existence and uniqueness of a mild solution provided, 
among other conditions, that 

fT 

/ \Z{t)\^^ dt < oo P-a.s. 
Jo 

The result could also be obtained applying theorem 10.15 of op. cit., provided one can 
prove that £ has cadlag trajectories in the domain of a fractional power of a certain 
operator defined in terms of A. In some specific cases, such condition is implied by 
suitable integrability conditions of the Levy measure. Unfortunately it seems to us rather 
difficult to verify such conditions, a task that we have not been able to accomplish. On 
the other hand, our approach, while perhaps less general, yields the well-posedness result 
under seemingly natural assumptions. 

Remark 3.9. By arguments similar to those used in the proof of theorem 13.61 one can 
also obtain that 

Esup|XA(t)P < C, 

t<T 

i.e. that {Xx} is bounded in L^(L°°([0, T]; 7^)). By means of Banach-Alaoglu's theo- 
rem, one can only conclude that X\ X in L^(L^([0, T]; 7^))', which is larger than 
L2(L°°([0,T];7i)). In fact, from ^ Thm. 8.20.3], being L^{[0,T];n) a separable Ba- 
nach space, one can only prove that if is a continuous linear form on L^(L^([0,T];7^)), 
then there exists a function / mapping 17 into L°°{[0,T];7i) that is weakly measurable 
and such that 

F{g)=E{f,g) 

for each g G ]L'^{L\[0,T];n)). 
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